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ABSTRACT 



Aims. The interpretation of hard X-ray emission from galaxy clusters is still ambiguous and different models proposed can be probed 
using various observational methods. Here we explore a new method based on Fe line observations. 

Methods. Spectral line emissivities have usually been calculated for a Maxwellian electron distribution. In this paper a generalized 
approach to calculate the iron line flux for a modified Maxwellian distribution is considered. 

Results. We have calculated the flux ratio of iron lines for the various possible populations of electrons that have been proposed 
to account for measurements of hard X-ray excess emission from the clusters A2199 and Coma. We found that the influence of 
the suprathermal electron population on the flux ratio is more prominent in low temperature clusters (as Abell 2199) than in high 
temperature clusters (as Coma). 

Key words. Galaxies: clusters: general; Atomic processes; Radiation mechanisms: non-thermal 
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1. Introduction 

Observations with BeppoSAX have detected hard X-ray tails in 
the X-ray spectra of some clusters as the Coma cluster (Fusco- 
Femiano et al. 1999) and Abell 2199 (Kaastra et al. 1998, 
Kaastra et al. 1999). These tails, which have been fit by power- 
law spectra, are in excess to the thermal bremsstrahlung X- 
ray emission from the hot intracluster medium (ICM). The ev- 
idence and the nature of hard tails in these and other clusters 
is discussed in the recent review by Rephaeli et al. (2008). 
The hard X-ray fluxes from galaxy clusters are usually inter- 
preted either as due to inverse Compton scattering (ICS) of rel- 
ativistic electrons on relic photons (Sarazin & Lieu 1998) or as 
bremsstrahlung emission from nonthermal subrelativistic elec- 
trons (see e.g. Sarazin & Kempner 2000) or from thermal elec- 
trons with a Maxwellian spectrum distorted by the particle ac- 
celeration mechanism (Dogiel 2000, Dogiel et al. 2007). 

The more traditional interpretation based on the ICS emis- 
sion from a relativistic electron population faces a serious prob- 
lem. The combination of hard X-ray and radio observations of 
the Coma cluster within the ICS model strongly indicates a very 
low magnetic field, B~0.1 fiG, much lower than the values de- 
rived from Faraday rotation measurements (see e.g. Clarke et al. 
2001). The situation in Abell 2199 is more extreme because no 
extended diffuse radio emission is detected from this cluster The 
discovery of the hard X-ray emission of the cluster A2199 im- 
plies a very weak ICM magnetic field of < 0.07yL(G if the hard 
X-ray emission is ICS (Kempner & Sarazin 2000). 

Bremsstrahlung radiation from suprathermal electrons with 
energies higher than 10 keV (nonthermal electrons or thermal 
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electrons with a distorted Maxwellian spectrum) may explain 
the hard X-ray excess emission observed in the Coma cluster 
and Abefl 2199 (e.g. Sarazin & Kempner 2000, Dogiel 2000). 
This is an alternative to the traditional but problematic inverse 
Compton scattering interpretation. Such subrelativistic electrons 
would form a particle population in excess over the thermal gas. 
A possible explanation for this population would be that they 
are particles being accelerated to higher energies, either by intr- 
acluster shocks or by turbulence in the ICM (e.g. Dogiel 2000). 

The best way to resolve the question of whether the observed 
hard X-rays are due to ICS or are evidence for a modified ther- 
mal distribution in clusters is to probe directly such a distribu- 
tion. 

The Sunyaev-Zel'dovich (SZ) effect signal, the spectrum of 
which depends on the electron distribution function in clusters 
of galaxies, can be used to discriminate among different inter- 
pretations of the X-ray excess (Colafrancesco 2007, Dogiel et 
al. 2007). The study of the influence of suprathermal electrons 
on the SZ effect was done for the Coma and Abell 2199 clusters 
by Blasi et al. (2000) and Shimon & Rephaeli (2002). However, 
realistical models of suprathermal electrons in Coma and Abell 
2199 predict a spectral distortion of the cosmic microwave back- 
ground radiation due to these electrons that is only a small frac- 
tion of the corresponding SZ effect due to the hot intracluster 
gas (see e.g. Shimon & Rephaeli 2002, Dogiel et al. 2007). 
Therefore the observation of the impact of suprathermal elec- 
trons on the cosmic microwave background will be challenging 
from the experimental side (see Dogiel et al. 2007). 

In this paper we consider a new probe to discriminate among 
different interpretations of the X-ray excess, namely the flux 
ratio of the emission lines due to FeKa transitions: FeXXV 
(helium-like) and FeXXVI (hydrogen-like). This flux ratio is 
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very sensitive to the population of electrons with energies higher 
than the ionization potential of a FeXXV ion (which is x 8.8 
keV) and is a promising tool to reveal the presence of suprather- 
mal electrons in galaxy clusters. 

A generalized approach to calculate the iron line fluxes is 
considered in Sect. 2. The flux ratio of the emission lines due 
to the FeKa transitions is calculated for the modified thermal 
distributions in the clusters Abell 2199 and Coma in Sect. 3. The 
possibility to separate the thermal and non-thermal components 
by using the shape of the bremsstrahlung continuum spectrum is 
discussed in Sect. 4. We draw our conclusions in the final Sect. 
5. 



2. The flux ratio of the FeXXV and XXVI iron lines. 

Since the fluxes of the FeXXV and FeXXVI Unes have the same 
dependence on the metal abundance, as well as on the emission 
measure, their ratio is independent of these parameters. This iron 
line ratio can be used to determine the temperature of the intra- 
cluster gas (e.g. Nevalainen et al. 2003). In this section we pro- 
pose a generalized approach to calculate the iron line flux ratio 
for modified Maxwellian electron distributions. 



2.1. Ionization and recombination rates. 

The ionization rates, recombination rates and emissivity in a 
spectral line have usually been calculated for a Maxwellian elec- 
tron distribution (e.g. Arnaud & Raymond 1992). However, in 
many low-density astrophysical plasmas, the electron distribu- 
tion may differ from a Maxwellian distribution. The influence of 
the shape of the electron distribution on the ionization and re- 
combination rates in various physical conditions was examined 
by Porquet et al. (2001). 

A Maxwellian distribution is generally considered for the 
electron distribution in galaxy clusters. Modified Maxwellian 
electron distributions expected in galaxy clusters with a hard X- 
ray excess seem to be reasonably described by a Maxwellian dis- 
tribution at low energy and by a power-law distribution at higher 
energy (e.g. Sarazin & Kempner 2000). 

It is convenient to express the electron distribution in term of 
the reduced energy x - E/kT: 



(1) 



where «e is the electronic density and k is the Boltzmann con- 
stant. 

Let us consider a collisional process with cross section o-{E), 
varying with the energy E of the incident electron. The corre- 
sponding rate coefficient F (cm^'s"') either for a Maxwellian dis- 
tribution or a modified thermal distribution, f{x), is obtained by 
averaging the product of the cross section by the electron veloc- 
ity over the electron distribution function: 



/2A:r^ 



\l/2 ^cx 



x^'^cr(xkT)f(x)dx 



(2) 



where nig is the electron mass, Xthi - Et^r/ikT), Ethr corresponds 
to the threshold energy of the considered process. 

For recombination processes, no threshold energy is in- 
volved and Xthi - 0. The rates are noted Ci, qrr and udr for the 
ionization, radiative and dielectronic recombination processes 
respectively, for a Maxwellian electron distribution. 

In equilibrium, the ionic fractions do not depend on the elec- 
tron density and the ionic fraction ratio ^peXxv/^FeXxvi of two 



adjacent stages FeXXV and FeXXVI for a Maxwellian distribu- 
tion can be expressed by: 



aR(FeXXVI) 
Ci(FeXXV) 



(3) 



where Ci(FeXXV) and Q'R(FeXXVI) are the ionization and total 
recombination rates of ions FeXXV and FeXXVI respectively. 

For the direct ionization cross section of FeXXV we use the 
parametric formula as in Arnaud & Rothenflug (1985): 



aui(E) = 



AjUj +BjUf + Cj In(Mj) -I- Dj ln(uj)/uj 



Mi/2 



(4) 



with u = E /Ij, Uj - I - 1/Mj, E is the incident electron energy, /j 
is the collisional ionization potential for the level j considered. 

The sum is performed over the subshells j of the ionized 
ion, for the ion FeXXV the Is subshell is considered (Arnaud & 
Raymond 1992). The parameters A, B, C, D (in units of 10 
cm^ eV^) and / (in eV) are taken from Arnaud & Raymond 
(1992). 

The autoionization process of the FeXXV ion can be ne- 
glected (Arnaud & Raymond 1992). 

The ratio of the ionization rate in a modified Maxwellian 
distribution over that in a Maxwellian distribution is: 



C x^^^o-ui{xkT)fMM(x)dx 
r°° x^'^o-Di(xkT)fM(x)dx 



(5) 



where /mm(^) is the modified Maxwellian distribution, fuix) is 
the Maxwellian distribution, ctdi is the direct ionization cross 
section of the ion FeXXV. The ionization rate is very sensitive 
to the fraction of electrons above the threshold energy. 

Recombination of a free electron can proceed either through 
a radiative free-bound transition (FeXXVI -H — > FeXXV + hv) 
or by a radiationless dielectronic recombination. 

The radiative recombination rates are less affected by a mod- 
ified thermal distribution than the ionization rates, since the cross 
section for recombination decreases with energy and there is 
no threshold. To estimate the radiative recombination rate ra- 
tio, we follow the method used by Owocki & Scudder (1983) 
and Porquet et al. (2001). The ratio of the radiative recombi- 
nation rate in a modified Maxwellian distribution over that in a 
Maxwellian distribution is: 



X ''fMM(x)dx 

Jo x-'TfM{x)dx 



(6) 



Following the method of Porquet et al. (2001) we use the value 
r] = 0.8 for an iron ion corresponding to the mean value < rj > 
reported in Arnaud & Rothenflug (1985). 

The dielectronic recombination is a resonant process involv- 
ing bound states at discrete energies £, and can be computed by 
summing the contribution of many such bound states. Following 
the method used by Owocki & Scudder (1983), we assume that 
the corresponding dielectronic recombination cross section can 
be approximated by: 



ctdr{E) = Yj^^^^^-^''^ 



(7) 



where D, are the dielectronic recombination coefficients. 
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The ratio of the dielectronic recombination rate in a modified 
Maxwellian distribution over that in a Maxwellian distribution 
is: 



j8dr = 



J^^ x^'^o-uR(xkT)fMM(x)dx 
x^'^o-BR{xkT)fM{x)dx 



(8) 



For the ion FeXXVI there is one bound state with the en- 
ergy El - 5.3 keV. If the break energy up to which the electron 
distribution is Maxwellian is higher than this bound state energy 
then the dielectronic recombination rate is not influenced by the 
suprathermal electrons. 

Thus the ionic fraction ratio ^peXxv/^FeXxvi of two adjacent 
stages FeXXV and FeXXVI for a modified Maxwellian distri- 
bution can be written as: 



ifreXXVI /mm 



/^RRcyRR(FeXXVI) +/jDRgDR(FeXXVI) 
y8iCi(FeXXV) 



(9) 



2.2. Excitation rates and iron line flux ratio. 

In the coronal model (see e.g. Mewe 1999), the line spectrum 
is dominated by radiative decay following electron impact ex- 
citation, plus a smaller contribution of recombination lines. We 
assume here that all iron ions that are to be excited are in the 
ground state (see e.g. Mewe & Gronenschild 1981). Considering 
only the dominant process of collisional excitation (Tatischeff" 
2003), the volume emissivity P^^^, (in units of photons cm"-' s"') 

of a particular line transition a ^ ^ in an ion Fe^' can be written 

as 



(10) 



where hh is the H ionic number density (cm^^*), ope is the abun- 
dance of iron relative to hydrogen, ^pg+, is the ionic fraction of 
ion Fe^', 5^"+, is the rate for electron impact excitation of an 

ion Fe+' from its ground state to its excited state a, and Bab is 
the radiative branching ratio of the transition a ^ b among all 
possible transitions from the level a. 

The excitation rates S p"^, {fix)} are functionals which are cal- 
culated by averaging the product of the corresponding cross sec- 
tion and electron velocity over electron distribution functions 
fix): 



5^';,[/(x)] 



™e / J.v„ 



x^l^alf(,xkT)f{x)dx , 



where cr^^ ' is the excitation section of the ion Fe^', Xthi- ex - 
Ethr,ex/{kT), and isthr.ex corresponds to the threshold energy of 
the excitation process. 

The emission lines due to FeKo- transitions of ions FeXXV 
and FeXXVI are at 6.7 keV and 6.9 keV respectively. Ions with 
closed-shell configurations are more stable than those with par- 
tially filled shells; thus He-like FeXXV, whose ground state is 
1 5^, is dominant in a large temperature range, because its ion- 
ization rate is relatively low compared to those of adjacent ions 
(e.g. Arnaud & Raymond 1992). The strongest line emission is 
then the He-like FeKa line complex at 6.7 keV that corresponds 
to transitions Is- - Islp^P, Lr - Islp^P, 1/ - ls2s^S. At 
high temperatures (e.g., kT - % keV), the hydrogen-like iron 
line (2p — > L? transition) at 6.9 keV also becomes intense. 

The electron impact excitation scaled cross sections for the 
helium-like ion FeXXV (including impact excitation from I'S' 



to 2'f , 2^P, 2^S levels) are taken from the article of Bazylev & 
Chibisov(1981): 



n ^ 19 im I 1-93 , 6.071n(z) \ , 
Z o-{ls -^ls2pP)-\ — - — h \naQ 



zV(li2 ^ ls2p^P) = ^^nal 

z 



zV(li^ ^ ls2s^S) = I — —\^aQ 



0.93 0.59 



(12) 
(13) 
(14) 



where aq is the Bohr radius, z is the incident electron energy in 
threshold units, Z is the ion nuclear charge. 

The electron impact excitation scaled cross section for the 
hydrogen-like ion FeXXVI is taken from the paper of Fisher et 
al. (1997): 



zV(ls 



,1.66 ln(z)\ 
2p) = (-5--H2.49— 1^ 



(15) 



There are also contributions from excitations to higher lev- 
els, that may radiatively decay to the upper levels of the He-like 
triplet and H-like doublet. These so-called cascade eff'ects gen- 
erally cannot be ignored. The electron impact excitation scaled 
cross sections of an iron ion from its ground state to the higher 
levels are taken from the article of Bazylev & Chibisov (1981). 

The flux ratio of iron lines FeXXV and FeXXVI taking into 
account electron impact excitation is then 



^ei 



P'-^(FeXXV)A4;^^^ 
pi-2(FeXXVI)A4;2^^j 



(16) 



where the volume emissivities P'"^(FeXXV) and 
P' "2 (FeXXVI) are for the He-like triplet and for the H- 
like doublet, and the energies AiSp^^^^y and AZSp^^^yj equal 6.7 
keV and 6.9 keV respectively. 

According to Eq. ( fTOl i the expression for the flux ratio in 
terms of ionic fractions and excitation rates writes as 



Re\ — 



"^FeXXV 



FeXXV 



Q 



1-2 

FeXXV 



A£, 



1-2 



^FeXXVI Q. 



1-2 

FeXXVI 



(17) 



where the rate coefficients are Gp^xxv 



' FeXXV 



(11) GpelxVI - 2 2 FeXXVI-^" 
a b{<a) 



a b(<a) 

The excited states b correspond 



to the upper levels of the He-like triplet and of the H-like dou- 
blet, and the radiative branching ratios are given by 



Bab - 



2 Aa 

c(<a) 



(18) 



All necessary transition probabilities Age are taken from 
Ralchenko et al. (2008). 

However, in plasmas in collisional ionization equilibrium, 
radiative recombination still contributes about 10% to the to- 
tal line flux. We calculated the rate coefficients for the contri- 
bution from radiative recombination to the spectral line forma- 
tion with equation (A. 9) in Mewe et al. (1985). The influence of 
the suprathermal electron population on the radiative recombi- 
nation rates is described by Owocki & Scudder (1983) (see Eq. 
|6]l. Although, as noted below, the ratio of the radiative recombi- 
nation rate /3rr depends slightly on the presence of suprathermal 
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electrons in the spectrum, the ionic fractions of FeXXVI and of 
FeXXVII and, therefore, the line emissivities change with the 
presence of suprathermal electrons. 

The line flux ratio taking into account both electron impact 
excitation and radiative recombination is given by 



^^^^^FeXXV ^ 
"^FeXXVI 



^1 



FeXXvGpeXXV 



+ freXXVlff] 



RR, FeXXV 



fexxviGpexxvi fexxviio;,' ^ 



(19) 



RR, FeXXVI' 



are the rate coefficients for the 



where p^^xv and a'^l p^^xvi 
contribution from radiative recombination to the spectral lines 
FeXXV (He-like triplet) and FeXXVI (H-like doublet) respec- 
tively. The ionic fractions of ^FeXxvii and ^FeXxvi were calculated 
following the same method as in Sect. 2.1. 

Figure 1 shows the iron line flux ratios for the pure 
Maxwellian distribution calculated in this section (dot-dashed 
line) and those obtained from the MEKAL model (points) by 
Nevalainen et al. (2003). 



3. An analysis of electron spectra in clusters. 

Bremsstrahlung from suprathermal electrons has been invoked 
as a possible explanation for hard X-ray tails in the X-ray spectra 
of some galaxy clusters. In this section we calculate the impact 
of suprathermal electrons on the FeXXV and FeXXVI emission 
line flux ratio in the clusters A2199 and Coma. 



3.2. The Coma cluster. 

The Coma cluster is a rich, hot, nearby (z=0.02) galaxy cluster 
Its average temperature is kT -d,.2 keV as derived from XMM- 
Newton observations (Arnaud et al. 2001). 

Hard X-ray radiation was detected in excess of thermal emis- 
sion in the Coma cluster by a first Beppo-SAX observation 
(Fusco-Femiano et al. 1999) and confirmed by a second inde- 
pendent observation with a time interval of about 3 yr (Fusco- 
Femiano et al. 2004). The reliability of the Fusco-Femiano et 
al. (1999, 2004) analyses was further discussed by Rossetti & 
Molendi (2004, 2007) and by Fusco-Femiano et al. (2007). 

The presence of a second component in the X-ray spectrum 
of the Coma cluster has also been derived from two RXTE ob- 
servations (Rephaeli & Gruber 2002). 

The spectrum of background and accelerated electrons was 
found by Gurevich (1960) from a kinetic equation describing 
stochastic particle acceleration: 



f'^^ (x) = 



2 



exp 



Jo 1 + ffg^ j I Jo 



jo 1 + ag^ 1^ 



(21) 



where the parameter cf = 9 ■ 10"^ was derived by Dogiel (2000) 
from the bremsstrahlung model for the origin of the hard X-ray 
emission from the Coma cluster. Here x - E/{kT) is the reduced 
energy and the integration has been done on the quantity g = 
p/ yImJiT that is the dimensionless momentum. 



3. 1 . The galaxy cluster Abell 2199 

A2199 is a bright cluster at redshift z=0.03. Its average gas 
temperature is kT - 4.7 keV (Kaastra et al. 1999). Spatially 
resolved spectroscopy shows a hard tail in the X-ray spec- 
trum of this galaxy cluster (Kaastra et al. 1998). To interpret 
this hard tail Sarazin & Kempner (2000) applied a non- thermal 
bremsstrahlung model with the electron distribution function 
/mmW given by: 

/mmW " /m(^), x<3 

/mmW = /m(^) + X > 3 (20) 

where yu = 3.33, X - 0.34 is found from the condition that the 
non-thermal electron population is 8.1% of the thermal popula- 
tion (Sarazin & Kempner 2000). 

The ratio of the radiative recombination rates jSrr (see Eq. 
|6ll is 1.013 for the electron distribution function /^'mW- Since 
the break energy of likT is higher than the bound state energy 
5.3 keV for the cluster temperature, the ratio of the dielectronic 
recombination rates is 1 . Therefore the recombination rates are 
not affected by suprathermal electrons. 

In Fig. 1 we compare the flux ratios R for a Maxwellian 
electron distribution (dot-dashed line) and for a modified 
Maxwellian distribution (dashed line) in the temperature 
range 4.5 keV < kT < 8.5 keV. 

For the galaxy cluster A2199 (kT=4.7 keV) the flux ratio 
R for a modified Maxwellian distribution /j^'^ decreases by » 
27% with respect to the case of a Maxwellian distribution. This 
value of the flux ratio would correspond to a thermal electron 
spectrum (i.e. without suprathermal electrons) with an effective 
temperature of kT = 5.4 keV. 
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T (keV) 

Fig. 1. Iron line flux ratios R for a Maxwellian electron distribu- 
tion (dot-dashed line), a modified Maxwellian distribution 

(dashed line) and a modified Maxwellian distribution f^^^ (solid 
line) in the temperature range 4.5 keV < kT < 8.5 keV 

The ratio of the radiative recombination rates ^Srr (see Eq.|6l) 

(21 

is 1.01 for the electron distribution function fy^y^x). For values 
kT in the range 4.5 - 8.5 keV, the values of ySoR are found in the 
range 1.0005-1.0025. Therefore the recombination rates are not 
affected by suprathermal electrons. 

In Fig. 1 we compare the flux ratios R for the Maxwellian 
electron distribution (dot-dashed line) and the modified 
Maxwellian distribution /j^j^ (solid line) in the temperature 
range 4.5 keV <kT < 8.5 keV. 

For the Coma cluster {kT-%.2 keV) the flux ratio R for a 
modified Maxwellian distribution /|^^ decreases by a; 9% with 
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respect to the case of a Maxwellian distribution. This value of the 
flux ratio would correspond to a thermal electron spectrum (i.e. 
without suprathermal electrons) with the effective temperature 
ofkT = 8.6keV. 

3.3. A synthetic low temperature cluster 

In Sects. 3.1 and 3.2 iron line flux ratios were calculated for low 
temperature and high temperature galaxy clusters (Abell 2199 
and Coma respectively). As shown in Fig.l, the impact of a 
suprathermal electron population on the iron line flux ratio is 
stronger in low temperature clusters. We demonstrate here for a 
specific example how the effective temperature inferred from the 
flux ratio of the iron lines can yield important constraints on the 
fraction of suprathermal electrons. For this purpose, a synthetic 
cluster with temperature kT - 4.7 keV and with an electron 
distribution function is considered. The dependence of the 
effective temperature on the fraction of suprathermal electrons is 
shown in Fig. 2. 
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-.5 ' ' ' ' ' ' ' ' ' 

5 10 15 20 25 30 35 40 
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Fig. 2. Dependence of the effective temperature kT (given in 
units of keV) on the fraction of suprathermal electrons in a clus- 
ter of temperature 4.7 keV. 

Since the effective temperature changes strongly with the 
suprathermal electron fraction, the iron flux ratio can be used 
to reveal a suprathermal electron population in low temperature 
clusters and in cluster cool cores. 

4. The continuum spectrum 

Non-thermal electrons change the flux ratio R and lead to an ap- 
parently higher temperature derived from the iron line ratio (i.e. 
the effective temperature, see Sect. 3 and Fig. 2). Alternatively, 
the temperature can be measured from the bremsstrahlung 
spectrum curvature. But, in the presence of non-thermal elec- 
trons, these same electrons will also give rise to non-thermal 
bremsstrahlung, therefore altering the shape of the continuum 
spectrum. A method to determine the temperature of the ther- 
mal part of a more complex electron distribution from the 
bremsstrahlung spectrum is proposed in this section. A disagree- 
ment between both temperatures will depend on the strength of 
the non-thermal electron component. 

To separate the contributions of thermal (low-energy) 
and of non-thermal (high-energy) electron components to the 



bremsstrahlung spectrum we study here the features of the en- 
ergy flux spectrum. 

The bremsstrahlung energy flux can be estimated as 

njV r°° [2E 
0(£,) = -^£, crB(E,E,)J—f(E)dE (22) 

where V is the cluster volume, d is the distance to the galaxy 
cluster, and o-BiE,E^) is the bremsstrahlung cross section. For 
the sake of illustration, let us consider the bremsstrahlung cross 
section in the form of cr^(E,E^ - const /E^E (following the 
approximation from the paper of Ginzburg 1979). Then, for a 
Maxwellian electron distribution, the dependence of the energy 
flux on the photon energy is given by 'J>{E^) cc exp{-E^/kT). 
For modified Maxwellian electron distributions, the energy flux 
at low photon energies has two components: thermal (i>i{E^) = 
Ci exp(-£'x/^7^) and non-thermal (^2{Ex) - C2. The non- 
thermal component of the bremsstrahlung energy flux is con- 
stant if the low limit of the integral in Eq.(l22li is smaller than 
the energy Em at which the electron distribution deviates from 
a Maxwellian distribution {Em = 14 keV for the cluster A2199 
(Kempner et al. 2000) and = 30 keV for the Coma cluster 
(Dogiel 2000)). Considering the energy band with E^ < Em, 
it is possible to fit the energy flux spectrum with a function 
(^{E^) - C\ exp(-E^/kT) + C2 and calculate the normalization 
constants Ci and C2, and also the temperature T of the thermal 
part of the electron distribution from the bremsstrahlung spec- 
trum. 

The analysis of RXTE measurements for the Coma cluster 
(Rephaeli & Gruber 2002) has shown evidence for the presence 
of a second spectral component at energies up to ~20 keV, since 
the fit to a single isothermal model has a poor quality. When a 
second thermal component is added, the best fit temperatures of 
the primary and secondary components are then kT[ - 7.5 keV 
and the very high value kT2 - 37.1 keV. The contribution of the 
second component (^2{Ex) exp(-£'x/^7^2) to the low-energy 
continuum spectrum is flat as noted above. Therefore the thermal 
and non-thermal components can, in principle, be separated by 
studying the shape of the continuum spectrum. 

From Eq.(l22li we calculated the energy fluxes <l>(£'x) of the 
Coma cluster in the range 4-20 keV for a Maxwellian electron 
distribution and for a modified Maxwellian distribution f^l^, as 
shown in Fig. 3. The values of the cluster parameters (tempera- 
ture, density etc) were taken from Dogiel (2000). The difference 
between the total 4-20 keV fluxes (i.e. the total 4-20 keV flux 
of the second spectral component) is x7%. The spectral compo- 
nents must be separated in order to obtain the temperature from 
the bremsstrahlung continuum spectrum. 

5. Conclusions 

We have shown in this paper that the iron line flux ratio de- 
pends on the presence of suprathermal electrons that have been 
proposed to account for measurements of hard X-ray excess 
emission from galaxy clusters. The influence of the energetic 
suprathermal electron population on the iron line flux ratio is 
more prominent in low temperature clusters (as Abell 2199) than 
in high temperature clusters (as Coma) because the fraction of 
thermal electrons with energies higher than the helium-like iron 
ionization potential in low temperature clusters is smaller than 
that in high temperature clusters. 

Since the decrease of the flux ratio of He-like to H- 
like Ka lines is expected for modified Maxwellian distributions 
in A2199 and Coma with respect to the case of a Maxwellian 
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Fig. 3. The energy flux <l'(£'x) of the Coma cluster in the range 4- 
20 keV for a MaxwelHan electron distribution (dashed line) and 
for a modified Maxwellian distribution (solid line). 



distribution, observation of the flux ratio is a tool to test the 
nonthermal electron bremsstrahlung model and to discriminate 
among different interpretations of the X-ray excess. In order to 
demonstrate the presence and to measure the strength of the non- 
thermal electron component we propose to compare the temper- 
atures obtained from the iron line flux ratio and from the low- 
energy continuum spectrum. 

The spectral resolution of XMM-Newton is sufficient to 
measure the flux ratio of the iron lines in hot temperature clus- 
ters. The constraint of the flux ratio for Coma within radius 5' 
is 1.6+[J-^ (Nevalainen et al. 2003). However, the XMM-Newton 
sensitivity in this high-temperature regime is insufficient to re- 
veal the contribution from suprathermal electrons in the Coma 
cluster (see Fig. 1). 

At low temperatures (e.g. kT < 5 keV) the FeXXVI line is 
weak and is within the noise level of the XMM-Newton data 
(Nevalainen et al. 2003). Therefore the iron line ratio cannot 
be measured by XMM-Newton in cooler clusters. On the other 
hand, the flux ratio of the iron lines in low temperature clusters 
can be detectable by XMM-Newton if the fraction of suprather- 
mal electrons is sufficiently high (see Fig. 2). 

Suzaku is also able to measure these two Fe lines in hot clus- 
ters due to its good spectral resolution (e.g. Fujita et al. 2008). 

We have considered the He-like triplet and the H-like dou- 
blet iron lines in this paper. Although XMM-Newton and Suzaku 
can distinguish the He-like from the H-like complex, their spec- 
tral resolution of ~100 eV causes that the observed line fea- 
tures do not only consist of the pure He-like triplet and H- 
like doublet, but each of these is blended with a multitude of 
satellite lines (e.g. Gabriel 1972; Dubau et al. 1981). For in- 
stance, for a temperature of 4.5 keV about ~30% of the flux 
from both line complexes is due to these satelUte Unes. In or- 
der to analyse the influence of satellite lines on the measure- 
ment precision of the iron line flux ratio we calculated the flux 
ratio Rb of the two (6.6-6.7 keV) and (6.9-7.0 keV) blends us- 
ing the line list which was taken from http://www.sron.nl/ divi- 
sions/hea/spex/versionl.lO/line/index.html and found that in the 
temperature range [4.5-8.5 keV] the values of R and Rb differ by 
less than 5%. Therefore the blend flux ratio Rb as well as the ra- 
tio R can be used to measure the temperature in this temperature 
range. 



The autoionizing levels responsible for the satellites are ex- 
cited by electrons at precisely the energies (see Eq.(40) from 
Mewe & Gronenschild (1981)) corresponding to those levels. 
Since the energies Eg are smaller than the energy at which 
the electron distributions in the Abell 2199 and Coma clusters 
deviates from a Maxwellian distribution, the exciting electrons 
belong to the thermal part of the electron distribution. Taking 
into account this fact and the dependence of the ionic fractions 
on the distribution function (see Sect.2.1) we estimate that the 
decrease of the flux ratio Rb for modified Maxwellian distribu- 
tions in A2199 and Coma with respect to a Maxwellian distribu- 
tion ai-e 30% and 10%. 

New high spectral resolution instruments with higher sensi- 
tivity such as XEUS are needed to resolve the lines and to mea- 
sure the flux ratio of the iron Ka lines for the purpose of testing 
the hard X-ray tail interpretations. 
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